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Abstract 

Differential cross sections for inclusive dijet photoproduction on a virtual pion have 
been calculated in next-to-leading order QCD as a function of Et,i], and Xt^. The 
cross sections are compared with recent ZEUS data on photoproduction of dijets 
with a leading neutron in the final state. 



1 Introduction 

Recently the ZEUS collaboration at HERA presented differential cross section data for 
the neutron tagged process e"^ + p — e"*" + n + jet + jet + X |l|]. The cross sections have 
been measured in the photoproduction region with photon virtuality < 4 GeV^ for jp 
center-of-mass energy W in the interval 134 < W < 269 GeV, for jets with transverse en- 
ergy Et > 6 GeV, neutron energy En > 400 GeV, and neutron production angle 6„ < 0.8 
mrad. The cross sections were compared to predictions with the one-pion exchange model 
and rather good agreement was found. 



Due to the kinematic constraints on the neutron detection the squared momentum trans- 
fer t between ingoing proton and outgoing neutron is very small [Q. In this case it is 
expected that the p ^ n transition amphtude, i.e. the amplitude for the production of 
neutrons, which carry most of the momentum of the proton, is dominated by the lightest 
particle in the t-channel, the pion. Hence the isolation of the one-pion exchange contri- 
bution provides the possibility to study the photon-pion interaction 7 + vr — A 0. 



Selecting high Et jets in the final state A allows us to apply the QCD improved parton 
model to describe the characteristic features of the jet production dynamics. This is com- 
pletely analogous to the process 7 + p — > jet + jet + X, which has been calculated in the 
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past up to next-to-leading order (NLO) using different computational techniques [Q, ^. 
The only difference in 7-|-7r jet-|-jet-|-X is that the parton distribution functions (pdfs) 
of the proton are replaced by the pdfs of the pion. In addition, the pion flux, generated 
by the p — > n transition, has to be supplied. 

The pion flux can, in principle, be measured in charge-exchange processes in soft hadronic 
reactions, where an initial-state proton is transformed into a final-state neutron, p — > n, 
with small momentum transfer. A successful phenomenological description of the corre- 
sponding data has been obtained in the framework of reggeized isovector exchange, such 
as TT, p, and 02, with the pion being by far the lightest meson dominating the p — >■ n tran- 
sition, particularly at small values of the squared momentum transfer t between proton 
and neutron p. 



The pdfs of the pion can be measured in deep-inelastic scattering with a forward going 
neutron detected in the final state: + p — > e'^ + n + X |^. The measurements of 
the structure functions so far are limited to rather small Bjorken-x 0, |^. The measured 
semi-inclusive cross sections for leading neutrons by the HI collaboration |§] can be de- 
scribed entirely by tt exchange, and the Bjorken-x and dependence of these data are 
consistent with leading order (LO) pdfs of the pion. There exist several such pion pdfs in 



the literature [^, |T0|, |TT| . They are constrained by dimuon and prompt-photon production 



data from fixed target experiments that are sensitive to the valence quark distribution in 
a Bjorken-x range relevant for dijet photoproduction on pions. 



In this work we present the results of a calculation of differential cross sections for the 
reaction e"*" + p — > e"*" + n -|- jet + jet + X with the kinematical constraints as in the ZEUS 
analysis We have done these calculations in LO and in NLO. It is well known, how- 
ever, that LO predictions are not reliable due to the possibly strong scale dependence of 
the results. Another problem of LO results for jet cross sections is the fact that in LO 
only two large jets with opposite transverse momenta can be produced. Therefore 
these cross sections are independent of the jet algorithm which is applied in order to de- 
fine jets in the analysis of the experimental data. In NLO calculations there appears an 
explicit dependence on the jet algorithm due to the possible production of three partons 
(in addition to the remnant partons on the photon or pion side, respectively) in the final 
state, two of which must be combined in one jet in accordance with the experimental 
constraints for the construction of jets out of hadrons. Therefore our LO results are just 
for comparison and for demonstrating the reduced scale dependence of the NLO results. 
Our genuine predictions are the NLO cross sections which will be compared with the 
ZEUS data ||. 

In the next section we shall describe the kinematical variables and define our input for 
the pion flux and the pion pdfs. Section 3 contains our results and the comparison with 
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the experimental data. In the last section we give a short summary and draw some con- 
clusions. 



2 Kinematical Variables and Further Input 

The event kinematics has been described in detail already in |T[]. Here we repeat those 
variables which are needed in the calculation of the cross cross sections. The reaction 

e+{k) +p{P) ^ e+{k') +n{P') + H, (1) 

where H is the hadronic system containing the jets, is characterized by the four-momenta 
k and k' of the initial and scattered positron and by P and P', the four-momenta of 
the ingoing proton and outgoing neutron, respectively. The positron-photon vertex is 
described by the exchanged photon virtuality and the positron's inelasticity y, i.e. 

= -q' = -{k-k')\ 

In the ZEUS experiment [[H < 4 GeV^ and 0.2 < y < 0.8, corresponding to the jp 
center-of-momentum energy range 134 < W < 269 GeV. The protons at HERA have 
the energy Ep = 820 GeV and collide with = 27.5 GeV positrons, corresponding to a 
center-of-momentum frame energy = 300 GeV. The two variables, which describe the 
proton-neutron vertex, are the fraction of the energy of the initial-state proton carried by 
the neutron xl and the square of the momentum transfer t between the proton and the 
produced neutron: 

P'k Er, 

XL ■ 



Pk Ep' 

t = f = {P- P'f. (3) 

The splitting function or pion flux for the transition p ^ + vr"*" is usually parameterized 
by different forms which can be summarized as 

UiA-L.t) = (4) 

Here gnnp is the coupling constant of the nirp vertex, is the pion mass, and a,r(^) = 
a'{t — m^) is the Regge trajectory of the pion. F{xL,t) is a form factor which describes 
the off-shell behavior of the virtual pion and/or possible final state rescattering effects 
of the neutron. Two choices for the form factor F{xL,t) used in the analysis of various 
charge-exchange scattering experiments are 

^ ^^^'^^ ~ 1 exp[R^{t - ml) /{I - xl)] ^ ' 
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where b and R are constants. The first of these choices, the exponential form, is usually 
taken with the Regge trajectory factor in (4) with a' = 1 GeV~^. The second choice, the 
light-cone form factor, is usually associated with the flux without the Regge trajectory 
factor, i.e. a' = in (4). These choices are associated with the experimental finding 
that in soft hadronic reactions the shape of the xl distribution depends on t |]T2[. The 
pion-nucleon coupling constant is well known from low energy ttN and A^A^ scattering 
data. We take 5'™/47r = 2 ■ 14.17 as obtained in a recent analysis [ITSj . 



In the scattering process of 2 — > 2 massless partons, the fractions of the four-momenta 
q = k — k' and q' = P — P' participating in the hard scattering by the initial state parton 
are given by 

where the sums in (6) and (7) run over the variables of the two jets in the final state. Here 
we assumed that q^ = q'^ = 0. The energy fraction contributing by the exchanged virtual 
photon to the production of the dijets is x^ whereas the corresponding contribution of 
the virtual pion (or possibly of a reggeized p or 02) is x^^. In (6) = yE^ is the energy 
of the ingoing virtual photon. In NLO also three jets can be produced in the final state. 
Then (6) and (7) are no longer valid. In order to estimate the energy fractions x^ and 
for this case, too, one uses (6) and (7) in the form that the sums in (6) and (7) run over 
the two jets of largest Et in an event. 



As is well known two mechanisms contribute to the photoproduction of jets, the direct 
and the resolved process. The observable x^ is sensitive to the amount of direct and re- 
solved processes. The LO direct process, where the photon couples directly to the quarks, 
contributes at x^ = 1, the resolved and the NLO direct processes contribute in the region 
x-y < 1. It must be emphasized, however, that the distinction of direct and resolved 
processes looses its meaning in NLO, so that the characterization of the contributions 
to the various x^ regions by the direct and resolved processes is not unique: It becomes 
scale dependent. The variable be used to analyze which region of the pion pdfs 

and which partonic contributions are most important for dijet production with leading 
neutrons. 



As for the jet definition and the combination of two partons into one jet in the NLO con- 
tributions of the direct and the resolved parts, we use the fcr cluster algorithm in the form 
as introduced in and as used in the jet analysis of the experimental data The 
photon flux is calculated with the usual Weizsacker- Williams approximation, including 
the non-logarithmic correction as calculated in |15|. For the resolved cross sections we 



4 



need the pdfs of the photon in LO and NLO. A popular parameterization is GRV |T6[. We 



use the parameterization GS96 ||T^, since it leads to a slightly better description of the 
experimental data. For the comparison with data of inclusive jet production without the 
leading neutron, which has been simultaneously analyzed in under similar kinematical 
conditions, we need also the pdfs of the proton, for which we employ the recent CTEQ5M 
||T8| parameterization. The A parameter which we need in the NLO a<j formula is taken 



from the proton pdf fit. It is A*^^^ = 326 MeV. For the pion pdfs we employ two alter- 



natives, SMRS3 and GRS ||TT|. Unfortunately, the A values used for the evolution of 
these pdfs are somewhat smaller, namely A*^^^ = 190 and 200 MeV. For LO predictions, we 
use LO matrix elements with the one-loop formula for as and the same values of A. Since 
only GRS have constructed LO and NLO pion pdfs, we have to use the NLO SMRS3 set 
also with our LO predictions. The GRS pdfs of the pion are constructed for three fiavors 
only. Since we include four fiavors throughout, we also need the pion pdf for charm. This 
is taken from the earlier work of GRV [DHIl. 



3 Results and Comparison with ZEUS Data 

The calculation of the cross sections is based on the formalism fully described in our pre- 
vious work 1^. To cancel infrared and coUinear singularities in the NLO contributions of 
the final state and to define the collinear initial state singularities to be removed through 
the renormalization of the pdfs of the photon or the pion, the phase-space slicing method 
was applied. Since the neutron kinematics could not be fixed in detail in the experiment 
we had to integrate over a finite region in and t. We did this in accordance with 
the specifications of the ZEUS experimental analysis. Except for the outgoing positron 
and the leading neutron the final state consists of two or three jets. The two-jet sample 
contains the bare parton jets from the LO and virtual NLO contributions and the two 
jets originating from the recombination of two partons in the three-parton terms. The 
three-jet sample is just given by the uncombined three-parton final states. 

First we calculated the differential cross section d'^a/dETdr], where Et and r] are the 
transverse energy and the rapidity of one of the jets in the two- or three-jet sample, re- 
spectively. As in the ZEUS analysis we restricted the selection in the three-jet sample. 
We included from this sample only the two jets with the largest Et- In our earlier work 
this cross section was denoted the inclusive one-jet cross section. There, however, we 
included all three jets from the three-jet contributions. In [|l| the cross section defined 
above is denoted the dijet cross section. This should be distinguished from the inclusive 
two-jet cross section defined in 0, which depends on Et, the transverse energy of one 
of the two jets with the highest Et, and the rapidities rji and ri2 of these high Et jets. 
The d'^a/dETdr] is integrated over — 2 < r/ < 2 as in [jl|. The result for the full jet sample 
(without detection of the leading neutron) as a function of Et for Et > 6 GeV turn out 
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to be in reasonable agreement with the experimental data from |I| (not shown). The 
point with the largest cross section at Et = 7 GeV agrees perfectly with the measured 
value. The points at medium Et lie somewhat higher than the theoretical curve. The 
equivalent cross section d?a / dETdrj, also integrated over |?7| < 2, for the case with the 
leading neutron is shown as a function of Et for 6.5 < Et < 29.5 GeV in Fig. la. In 
this figure three curves are shown: one from LO with pion pdfs SMRS3 and photon 
pdfs GS96 ll^, the latter in LO, and two other curves for NLO with pion pdfs SMRS3 



and GRS |TT] and charm content from GRV [10| for comparison. We see that all three 



predictions lie very near to each other in the logarithmic plot. The LO and NLO curves 
for SMRS3/GS96 differ only very little in this plot, showing that the NLO corrections are 
small for the scale choice /i = max(£'7^ i, Et,2) as used in Fig. la (and in all the following 
figures) for the LO and NLO result. The GRS choice for the pion pdfs leads to a some- 
what larger cross section at small E't's and smaller cross section at the larger Et^s as 
compared to the SMRS3 choice. We have calculated that the gluon in the pion gives the 
dominant contribution to the cross section for Et < 22 GeV. Resolved photons are more 
important than direct photons for Et < 9 GeV. The experimental data [|T| with statistical 
and systematic errors shown separately come with a light shaded band which indicates 
the additional systematic experimental uncertainty due to the calorimeter energy scale 
|l|]. The medium shaded band gives the uncertainty of the LO result when the scale of 
as and of the factorization is varied in the interval between Et/^ and 2 • Et- The scale 
uncertainty of the NLO prediction (dark shaded band) coincides with the line thickness 
of the NLO curve. Taking into account the fairly large experimental errors including the 
calorimeter energy scale uncertainty the agreement of the ZEUS data with the theoretical 
predictions concerning shape and absolute normalization is reasonably good. It must be 
emphasized, however, that the absolute normalization depends on the pion form factor 
(5). In Fig. la we assumed the light-cone form factor with R = 0.5 GeV^^. Recent 
determinations of from a comparison to the rather accurate neutron production data 
from Flauger and Mdnning yield R^ = 0.2 GeV-^ [|| and R^ = 0.4 GeV'^ [Tg] in rea- 
sonable agreement with our value for R^. A change of R from R = 0.5 to R = 0.6 GeV~^, 
i.e. an increase by 20%, leads to a 15% smaller cross section. 



Since we have calculated also (Pa/dETdrj with the same kinematical constraints for 
e + p — > e' + 2 jets + X we can take the ratio of the two cross sections da/dET for 
e + p^e' + n + 2 jets + X and e + p ^ e' + 2 jets + X. This is shown in Fig. lb and 
compared to the ZEUS data for this ratio. The agreement with the data is perfect. For 
low Et, where the experimental errors are smallest, the SMRS3 result agrees better than 
the ratio obtained with the GRS pion pdfs, otherwise the results are the same. In Fig. lb 
the experimental uncertainty from the calorimeter energy measurements, the systematic 
error and the scale uncertainty in the theoretical curve are not included since they are 
supposed to cancel in the ratio |]l|. Thus, this ratio seems to be a good example to get 
information on the pion pdfs and/or the pion form factor. 
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ep ^. e' + n + 2 jets + X at Vs = 300 GeV 
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Figure 1: Dependence of the dijet photoproduction cross section with a leading neutron 
on the transverse energy of one of the two jets (a). Two NLO predictions with different 
pion structure functions and one LO prediction are shown together with experimental 
data from ZEUS with statistical (inner) and systematic (outer) experimental error bars. 
In addition, error bands from the theoretical scale uncertainty (medium: LO, dark: NLO) 
and the experimental energy scale uncertainty (light) are shown. The NLO error band 
coincides with the thickness of the full line. In b) we present the ratio of the leading 
neutron over the inclusive cross section versus the transverse energy. NLO calculations 
with two different pion structure functions are shown together with the experimental 
ZEUS data. In this case, only statistical error bars are given. 
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The next calculated differential cross section is da/dt], i.e. d'^a/dETdrj integrated over 
Et with -ET,min = 6 GcV. The result is shown in Fig. 2a, where da/drj is plotted as 
a function of rj for 7] between —2 and 2, together with the data points from There 
are three theoretical curves for the three cases as in Fig. la. Here one can see the 
difference between SMRS3/GS96 in LO and NLO. The difference is very small, the LO 
cross section is approximately 5% larger than the NLO cross section for the same scale 
/i = max(i?r 1, i?T 2)- Concerning the dependence on the pion pdfs the choice SMRS3 
agrees best with the data. The result with the GRS pdfs is larger, in particular for 
negative 77. In two other pion pdfs have been published, SMRSl and SMRS2. The 
cross sections for these two pdfs lie between the SMRS3 and the GRS curves and are 
therefore not shown explicitly. The gluons in the pion contribute about twice as much 
as quarks to the rapidity distribution over the entire ?7-range. Direct (resolved) photons 
dominate in the region < (77 > 0). All three curves do not agree well with the data 
points for rj < —0.5. For jet production without leading neutron this region is known to 



be influenced by hadronization corrections [20|. We expect such corrections also for the 



cross section in Fig. 2a, although somewhat less in the case with leading neutron due 
to the reduced total energy producing the hard jets. The theoretical result with such 
hadronic corrections included (taken from the case without leading neutron |^) is shown 
in Fig. 2a as the lower full curve. It leads to a better agreement in the negative rj region, 
but it produces also a somewhat lower cross section at positive r^'s, reducing the good 
agreement wit the data in this region. In Fig. 2a we show again three shaded bands: (i) 
the light shaded band shows the additional systematic uncertainty from the calorimeter 
energy scale, (ii) the medium shaded band gives the variation of the LO result (dotted 
curve) with the scale fi, which was varied as usual in the region Ex/2 < fi < 2 ■ Et, 
(iii) the dark shaded band shows the variation of the NLO cross section with changing 
the scale in the same interval. Here, we see quite clearly that the NLO result has a very 
much reduced scale dependence compared to the LO result. Similarly to Fig. lb we have 
plotted in Fig. 2b the ratio of the cross sections for SMRS3/GS96 and GRS/GS96 and the 
cross section da/drj for dijet production without leading neutron as a function of rj and 
compared them with the same experimental ratio from ZEUS. Within errors the SMRS3 
and the GRS results are consistent with the data, the SMRS3 curve agrees slightly better. 



The last differential cross section which we investigated is da/d\ogiQ{x.,r) as a function 
of lo gio(^7r)- In order to determine Xtt, which is a measure of the exchanged pion's mo- 
mentum fraction participating in the production of jets as defined in (7), we need the 
transverse momenta and rapidities of the two most energetic jets in the event. The ex- 
perimental cross section da/dlogiQ^XT,), however, has the constraint that for the two jets 
Et,i > 6 GeV and i?^- 2 ^ 6 GeV. Cross sections for E^.i = -E't,2 are problematic theoret- 
ically, since they become infrared sensitive in NLO, i.e. depend on the slicing cut used 
to cancel infrared and collinear singularities. In order to avoid this sensitivity and also 
possible hadronic and intrinsic /ct effects not considered in our NLO framework we need 
constraints on Et,i,Et^2 or Et^s which smear out the problematic region Et,i = -Et,2- 
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ep ^. e' + n + 2 jets + X at Vs = 300 GeV 
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Figure 2: Dependence of the dijet photoproduction cross section with a leading neutron 
on the rapidity of one of the two jets (a). In b) we show the ratio of the leading neutron 
over the inclusive cross section versus the rapidity. Details are as in Fig. 1. The lower full 
curve in a) demonstrates the influence of hadronization corrections. 
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Figure 3: Dependence of the dijet photoproduction cross section with a leading neutron 
on the logarithm of the observed momentum fraction of the partons in the pion log^gl^Tr)- 
Details are as in Figs, la and 2a. In addition, we show the effect of lowering the transverse 
energy difference A from 1 GeV to 0.5 GeV. 



This problem is well known and was encountered some time ago also in connection with 



the calculation of the inclusive dijet cross section in 7j9 collisions One remedy to 

remove the infrared sensitivity is the requirement of slightly different lower limits on Et^i 
and Et,2- Therefore, we have calculated da/dlogiQ^x-jr) with the constraint Et,i > 6 GeV, 
Et,2 > 5 GeV if Eta > Et,2 or Et^i and Et,2 interchanged if Et,2 > -Er,i- The cross 
sections calculated with this constraint are plotted in Fig. 3, again for the three cases as 
in Fig. la. These cross sections are compared with the corresponding ZEUS data. The 
NLO prediction with the pion pdfs of SMRS3 agrees better with the data than the one 
with GRS pion pdfs. Except for large > 0.35 the gluons in the pion dominate over 
quarks, and the resolved photon contribution is larger than the direct one for > 0.04. 
The agreement with the data is good for x,r > 0.06, but less good for the smaller x-,,- This 
region is related to the negative t] region in the rj distribution in Fig. 2a, where one has 
hadronization corrections (see (7) for the relation between x,r and the ?7j's). We note that 
the K factor between the NLO and LO result is larger than in da/drj and is larger than 
one now in particular for small x^. This is connected with the specific constraint on Et,i 
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and Et,2, which was not necessary for the results in Figs. 1 and 2. This constraint does 
not exactly correspond to the one applied to the measured da/dlogiQ^XT,). On the other 
hand the low Et calorimeter jets used in the experimental analysis have an uncertainty of 
5% which introduces a corresponding uncertainty on the constraint Et,i,Et,2 > 6 GeV. 
In order to see the influence of the Eta,Et,2 cut we reduced the difference A between 
and i?T,2 from 1 GeV to 0.5 GeV. The result is shown in Fig. 3 as the lower full 
curve. The reduction of the difference leads to a lower cross section as we expect, so that 
the low experimental points are described very well now but not the higher x^^ region. 
Since this region is much less affected by hadronization corrections we prefer the original 
prediction (upper full curve in Fig. 3). The absolute normalization of the cross section 
da/dlog^Q^Xj^) depends also on the pion form factor. The curves in Fig. 3 are calculated 
with the light-cone form factor with R = 0.5 GeV^^. R = 0.6 GeV~^ instead gives a 15% 
lower cross section and better agreement for x^r < 0.06 but worse agreement above this 
value. Again taking into account possible hadronization corrections we consider the value 
R = 0.5 GeV~^ as more realistic. The result with R = 0.6 GeV~^ is almost indistinguish- 
able from the result with the exponential form factor with 6 = 0. 

In several uncorrected distributions are presented, in particular event rates as a function 
of x^ in order to say something about the amount of resolved and direct photoproduction. 
From these rates, the ratio of resolved (corresponds to x^ < 0.75) to the direct cross 
section (corresponds to x^ > 0.75) has been deduced. We checked that the theoretical 
prediction for this ratio is in reasonable agreement with the data. It decreases with in- 
creasing Et, as one expects it. 

In this work we considered only the pion exchange in the t-channel. Other isovector 
exchanges, such as p and 02, are expected to be small. Recent studies of leading neutron 
production in ep collisions show that the rate can be described entirely by 7r+ exchange 
0. The additional exchanges, which increase the rate of neutron production, are offset 
by absorptive rescattering of the neutron, which decreases the rate by approximately 
the same amount 0]. Therefore both effects are neglected in the present analysis. If 
the experimental data were more specific concerning the dependence on t and xl, the 
dominance of the vr"*" exchange could be demonstrated. 

4 Summary and Concluding Remarks 

We have presented a next-to-leading order calculation of dijet production in collisions of 
almost real photons and pions, where the photons are generated through bremsstrahlung 
from an initial electron or positron beam and the pion is assumed to give the dominant 
contribution to the forward proton to neutron transition amplitude. Our results have 
been compared with recent data from ZEUS for distributions in the transverse energy 
Et and rapidity r] of one of the two jets and for the observed momentum fraction of the 
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partons in the pion x-^. We found that the NLO corrections to the LO result were gen- 
erally small, also in the case of inclusive photon-proton scattering with similar kinematic 
constraints, but lead to a largely reduced dependence on the renormalization and factor- 
ization scales. For the NLO distribution in logio(a^7r), we had to relax the experimental 
condition of Et,i, -Et,2 > 6 GeV to allow for a small difference A = 0.5...1 GeV of the two 
minimal transverse energies, corresponding to the non-zero experimental uncertainty in 
the jet energy scale. 

The overall normalization of all three distributions is sensitive to the form factor of the 
p n7r~^ transition. We have used a light-cone form with a value for the parameter R 
which agrees well with recent determinations from neutron production data. Once the 
normalization is fixed, the shapes of the distributions, particularly the one in logiQ^x-,,), 
are sensitive to the parton distributions in the pion. The ZEUS data for the distribution in 
log;^Q(x^) as well as in Et and rj are best described by the SMRS3 pion parameterization, 
whereas the GRS parameterization with GRV charm distribution gives larger predictions, 
especially at small Xtt. We emphasize that for small Et one is particularly sensitive to 
the gluon content of the pion. Furthermore, the GS96 photon densities agree better with 
the data than GRV. Hadronization corrections were found to reduce the perturbative 
prediction considerably in the region of small 1], which corresponds to small Xt^. 



Acknowledgments 

We thank M. Khakzad, G. Levman and J. Whitmore for making the ZEUS data available 
to us and M. Wing for calculating the hadronization corrections with the HERWIG Monte 
Carlo generator. Financial support by the Deutsche Forschungsgemeinschaft through 
Grant No. KL 1266/1-1, by the Bundesministerium fiir Bildung und Forschung through 
Grant No. 05 HT9GUA 3, and by the European Commission through the Research Train- 
ing Network Quantum Chromodynamics and the Deep Structure of Elementary Particles 
under Contract No. ERBFMRX-CT98-0194 is gratefully acknowledged. 



References 

[1] J. Breitweg et al. [ZEUS Collaboration], DESY-00-142, October 2000, ^ep 



ex/OOlOOlS . 



[2] G. Kramer, Proc. of the Seminar on e-p and e-e Storage Rings, ed. by J.K. Bienlein, 
I. Dammann and H. Wiedemann, DESY-73-66, December 1973. 

[3] M. Klasen and G. Kramer, Z. Phys. C 76 (1997) 67; 

M. Klasen, T. Kleinwort and G. Kramer, Eur. Phys. J. direct C 1 (1998) 1. 



12 



[4] B. W. Harris and J. F. Owens, Phys. Rev. D 56 (1997) 4007; 
S. Frixione and G. Ridolfi, Nucl. Phys. B 507 (1997) 315; 

P. Aurenche, L. Bourhis, M. Fontannaz and J. P. Guillet, Eur. Phys. J. C 17 (2000) 
413. 

[5] U. D'Alesio and H. J. Pirner, Eur. Phys. J. A 7 (2000) 109 and the earher hterature 
given there. 

[6] J. D. Sulhvan, Phys. Rev. D 5 (1972) 1732; 

N. S. Craigie and G. Schierholz, Nud. Phys. B 100 (1975) 125. 

[7] M. Derrick et al. [ZEUS Collaboration], Phys. Lett. B 384 (1996) 388. 

C. Adloff et al. [HI Collaboration], Eur. Phys. J. C 6 (1999) 587. 

[9] P. J. Sutton, A. D. Martin, R. G. Roberts and W. J. Stirhng, Phys. Rev. D 45 (1992) 
2349. 

[10] J. F. Owens, Phys. Rev. D 30 (1984) 943; 

P. Aurenche, R. Baier, M. Fontannaz, M. N. Kienzle-Focacci and M. Werlen, Phys. 
Lett. B 233 (1989) 517; 

M. Gliick, E. Reya and A. Vogt, Z. Phys. C 53 (1992) 651. 

[11] M. Gliick, E. Reya and I. Schienbein, Eur. Phys. J. C 10 (1999) 313. 

[12] J. Engler et al, Nucl. Phys. B 84 (1975) 70; 

W. Flauger and F. Monnig, Nucl. Phys. B 109 (1976) 347; 
B. Robinson et al, Phys. Rev. Lett. 34 (1975) 1475; 
G. Levman and K. Furutani, DESY-95-142, July 1995, 
and further references given there. 



[13] T. E. Ericson, B. Loiseau and A. W. Thomas, |hep-ph/0009^ 



[14] S. Catani, Y. L. Dokshitzer, M. H. Seymour and B. R. Webber, Nucl. Phys. B 406 
(1993) 187; 

S. D. Ellis and D. E. Soper, Phys. Rev. D 48 (1993) 3160. 
[15] S. Frixione, M. L. Mangano, P. Nason and G. Ridolfi, Phys. Lett. B 319 (1993) 339. 
[16] M. Gliick, E. Reya and A. Vogt, Phys. Rev. D 46 (1992) 1973. 
[17] L. E. Gordon and J. K. Storrow, Nucl. Phys. B 489 (1997) 405. 
[18] H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 12 (2000) 375. 
[19] H. Holtmann, A. Szczurek and J. Speth, Nucl. Phys. A 596 (1996) 631. 



13 



[20] B. W. Harris, M. Klasen and J. Vossebeld, Proc. of the Workshop on Monte Carlo 
Generators at HERA Physics, Hamburg, Germany, 1998, |hep-ph / 9905348| ; 
M. Wing, priv. communication. 

[21] M. Klasen and G. Kramer, Phys. Lett. B 366 (1996) 385. 



14 



